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Formulae  are  given  for  the  near  field  pressures,  far  field 
pressures  and  acoustic  particle  displacements  of  a  ring  source 
and  a  ring  axial  force.  These  quantities  are  necessary  ingredients 
for  calculations  of  far  field  radiation,  in  the  presence  of  the 
scattering  action  of  an  axi symmetric  shell  structure.  Numerical 
examples  are  given.  _________ _ _ ^  j 
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1.  INTRODUCTION 


Previous  work  C13  has  presented  a  formulation  for 
calculating  numerically  the  far  field  scattering  hy  a  finite 
elastic  axisymmetric  shell  structure  which  is  subject  to 
general  external  acoustical  excitation,  field  quantities  being 
represented  by  Fourier  series  expansions  in  the  circumferential 
coordinate.  Specifically,  the  work  shows  that  if  expressions 
are  available  for  the  pressure,  particle  displacements  and  far 
field  pressure  of  the  excitation,  in  the  absence  of  the  shell, 
then  the  necessary  ingredients  are  present  for  numerical 
calculations  of  the  far  field  pressure  in  the  presence  of  the 
shell  Time-harmonic  plane  wave  excitation,  the  basis  of  most 
scattering  studies,  demonstrated  the  usefulness  of  the  method. 

-NIn  many  hydroacoustics  problems  the  excitation  is  located 
close  to  the  shell's  surface:  it  may  be  of  monopole  type  due  to 
cavitation  sources,  dipole  type  due  to  fluctuating  fluid  forces 
or  quadrupole  type  due  to  turbulence  sources.  Herein  the  excit¬ 
ation  cases  of  an  axial  ring  source  and  an  axial  ring  force  are 
analysed.  They  model,  in  a  very  simple  way,  propulsor  generated 
cavitation  and  fluctuating  thrust  in  the  surrounding  fluid.  The 
presence  of  a  shell  is  important  because  interaction  between  the 
source/force  and  a  resonant  scatterer  may  result  in  con&itferable 
amplification  or  attenuation  of  the  sound  field  in  certain  spa¬ 
tial  directions. 


2.  RING  SOURCE /FORCE  EXCITATION 


(a)  General 


Figure  1  shows  an  axisymmetric  elastic  shell  which  is 
excited  by  a  time-harmonic  ring  source  or  axial  force.  In  order 
to  predict  the  far  field  sound  radiation  using  a  coupled  finite 
element  and  Helmholtz  integral  equation  formulation  described 
elsewhere  C13  it  is  first  necessary  to  calculate,  in  the  absence 
of  the  shell,  the  near  field  pressures,  the  near  field  particle 
displacements  and  the  far  field  pressures. 


(b)  Free-Soace  Pressure 

The  free-field  pressures  of  an  off-axis  unit  point  source 
and  axial  force,  in  a  fluid  of  sound  velocity  c  and  density  p, 
are  well  known,  viz. 

ps(r,$,z)  =  exp( ikD0 ) /DQ 

(2.1) 

p^(r,$,z)  =  -  ( l/4ir )  (  3  /  3z  )  Cexp(  ikDp ) /Dq] 
where 

D2  =  r2+rQ+(z-z0)2-2rr0cos($-#0)  (2.2) 

iz  the  square  of  Lh;  distance  between  the  source/force  located 
at  (r0,$0,zQ)  and  the  observation  point  located  at  (r,$,z):  'k' 
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Is  the  acoustic  wavenumber,  w/c.  Time  variation  exp(-iwt)  is 
omitted  throughout. 

Expand  these  pressures  as  Fourier  series  in  the  circum¬ 
ferential  coordinate  difference  <P~4>0,  viz. 

Ps(r,$,z)  =  E~=0  psn(r,z)cosCn(iti-(ti0)3 

(2.3) 


pf(r,<|),z)  =  e"_0  pfn(r,z)cosCn(i(i-iJ)0)3 


where 

pgn(r,z)  =  (en/2ir)  /q*  Cexp(  ikDg ) /DglcosCnf  4>-<()0 )  3d<)> 

(2.4) 

pfn*r,z*  =  -  ( l/4ir )  ( en/2ir ) (  3  /  3z )  Cexp(  ikDg )  /Dg]cosCn(  <J> — )  3d$ 

in  which  e^=l  if  n=0  and  en=2  otherwise.  Using  the  periodic 

property  of  the  cosine  function  it  is  not  difficult  to  show  that 
these  equations  are  independent  of  $0.  Thus 

Psn<r,z)  =  (en/ir)f  *  Cexp(  ikD) /DDcosfntyldiJ) 

(2.5) 

pfn(r'z)  =  (en/lrl /q  (  3/8z>  Cexp<  ikD>  /D3cos(ntJ))diJ) 

where 

=  r^+r^+t z-z0 > 2-2rrQCOS ( xp  >  (2.6) 


Consider  a  ring  source/force  of  radius  rQ  and  amplitude 
B( <t>Q )  / 2irr 0  per  unit  length  of  the  circumference:  the  pressure 
fields  may  be  synthesized  as 

P8<r,0,z)  =  l“=0  anpgn(r,z)cos(n0) 

+  En=0  bnpsn(r'z)sin!n$) 

(2.7) 


pf  (  r ,  $ ,  Z  ) 


=  En=0  anPfn(r,z)cos(n$) 

+  En=0  Vfn(r'z)sin<n<l) 


where 


4 


(2.8) 


an  =  /0  CB($0>  /2irr0Dcos(n$ri)rrid$ 


0'‘0  *0 


bn  =  /q  CB($q) /2irr0Dsin(n(t>0)r0d(ji0 


Specializing  to  the  case  of  constant  force/source  amplitude,  per 
unit  length,  gives  the  ax i symmetric  contribution  from  the  n=0 
harmonic  alone,  viz. 

Ps(r,z)  =  Bd/irJJ'Q  Cexp<  ikD)  /D3dty 

(2.9) 

pf(r,z)  =  ( -B/4ir)  ( 1/ir) (  3/3z  )  Eexpl  ikD) /D3dip 

As  a  consistancy  check  let  rp-+0  in  equations  (2.1)  and  (2.9)  to 

give  identical  results  for  the  pressure  of  a  point  force/source 
on  the  axis  of  symmetry. 


< c )  Derivatives 

The  derivatives  required  in  the  analysis  are  obtained  by 
application  of  the  formulae 

3/3z  =  (3D/3Z) (d/dD)  =  C(z-z0)/DJ(d/dD> 


3/3r  =  (  3D/3r)  (d/dD)  =  C < r-rQcos ( i|>)  )  /D3(d/dD) 


(2.10) 


(d/dD)Cexp( ikD) /D]  =  Cik/D-l/D  3exp(ikD) 


giving 


(3/3z)Cexp(ikD)/D3  =  ( z-z0 > ( ik/D2-l /D3 >exp( ikD) 

( 3/3r ) Cexp( ikD) /DJ  =  C(  r-rpCos  ( i|i) )  3(  ik/D2-l/D3 >exp(  ikD) 


(2.11) 


(32/3z2)Cexp(ikD)/DD  = 


( ik/D2-l /D3 )exp( ikD)  +  ( z-zQ ) 2 ( -3ik/D4+3/D5-k2/D3 ) exp( ikD) 


<  3  /3r3z ) Cexp( ikD) /D]  = 

(z-zn)  (r-r.cosOli)  ) (-3ik/D3+3/D4-k2/D2)exp< ikD) 


(d)  Acoustic  Particle  Displacements 

The  acoustic  particle  displacement,  in  the  direction  normal 
to  the  elastic  surface,  of  a  pressure  field  p  is 
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(2.12) 


W(r,$,z)  =  Ccos((3>3p/3r  -  sin(0) 3p/3z3/pw2 

where  g  is  the  normal  angle  shown  in  Figure  1.  The  displacements 
for  the  source  and  force  excitation  are  found  by  replacing  'p' 
by  ps  and  p^  of  equations  (2.7). 

For  the  free  field  ring  source  the  required  derivatives  of 
the  Fourier  amplitudes  are 

3psn/3r  =  (en/ir >/*  (  3/3r )  Cexp(  ikD) /D3cos  (n\J>)dty 


(2.13) 

3psn/3z  =  (en/ir)/Q  (  3/3z )  Cexp(  ikD) /D3cos  (ntj>)dty 
and  for  the  free  field  ring  force  they  are 

3Pfn/3r  =  ( -l/4ir )  ( en/ir ) (  32/3z3r  )  Cexp(  ikD) /D3cos  (nij>)dt|) 

(2.14) 

3pfn/3z  =  (-l/4ir)  (en/ir )/*  (  32  /  3z2  )Cexp(  ikD'  /D3cos  ( mp)  dip 
Equations  (2.11)  give  derivatives  of  exp(ikD)/D. 


(e)  Far  Field  Pressure 

The  distance  factor  'D‘ ,  given  by  equation  (2.6),  is 
written  as 

D  =  (  r2+z2+rQ+ZQ-2zz0-2rr0cos(i|)>  ) 1 /2 

=  RCl+(rQ+ZQ)/R2-2zz0/R2-2rr0cos(iJ>>/R231/2  (2.15) 

In  the  acoustic  far  field,  R-*»,  the  second  term  is  neglected, 
giving  to  a  first  order  of  accuracy 

D  =  R(l-zocos(9)-rosin(0)cos(il>) )  (2.16) 

Substitute  equation  (2.16)  into  equations  (2.5)  to  give 
the  approximations, 

Pan(R,0>  =  (en/jr)Cexp(ikR) /RDexpt -ikzocos(0) )  X 
exp(-ikrosin(0)cos(ij>)  )cos(nty)d\l> 

(2.17) 

pfn(R,0)  =  <en/*)Cexp<ikR>/R]exp<-ikzocos<0>'  X 

(-ik/4ir)cos(0) /q  exp<-ikrosin(0)cos(t(») ) cos  ( nt|>)d\|> 

The  integral  can  be  evaluated  as  it  is  a  standard  definition  of 
the  Bessel  function 
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(2. IB) 


Jn<x)  =  (i_n/ir)/*  exp(ix.  cos(iJ>)  )cos(nT(i>di|> 

Thus 

pgn(R,8)  =  enCexp( ikR) /R3exp( -iks^cos ( 8 ) ) inJn( kr^sinCQ ) ) 

(2.19) 

pfn(R,6)  =  (  -lk/4if )  cos ( 6  )  X 

e^Cexpf  ikR)  /R3exp(  kZgCos ( 6)  )  inJn( kr^sinO )  ) 

The  far  field  pressures  for  general  amplitude  variation  of  the 
excitation  are  then 

P  (R,0,<t>)  =  £°°  n  a  p  (R,8)cos(n$! 

3  n=u  n  sn 

+  T°°  n  b  p  (R,6)sin(n$) 
n=0  n  sn 

(2.20) 

pf(R ,©,<!>)  =  l"_0  anpfn(R,8)cos(niti) 


+  rn=0  bnpfn(R'8)sin(n<,l) 

where  an  and  bn  are  again  defined  by  equations  (2.8).  Far  field 

pressures  for  constant  amplitude  of  the  excitation  are  those  of 
the  n=0  harmonic  alone,  viz. 

Ps (R,0)  =  BCexp( ikR) /R3exp(-ikzQcos\  Jo(krosin(0) ) 

(2.21) 

pf (R ,0 )  =  <-ikB/4ir)cos(8)  X 

Cexp( ikR ) /R3exp( -ikz0cos(8) ) J0(kr0sin( 0) ) 

In  the  low  frequency  regime  kr0<<l  the  Bessel  function  may 
be  replaced  by  its  value  for  a  small  argument.  In  this  case  the 
far  field  pressure  is  proportional  to  (kr0>n  which  is  independent 

of  the  ring  radius  when  n=0,  and  drops  rapidly  as  'n'  increases. 
However,  the  near  field  pressure  and  hence  the  excitation  are 
usually  strongly  dependent  on  the  radius. 


3.  NUMERICAL  TESTS 

A  Fortran  program  C13  for  calculating  far  field  scattered 
pressure  when  the  excitation  is  a  plane  wave  at  an  arbitrary 
incidence  has  been  modified  to  calculate  far  field  pressure  when 
the  excitation  is  a  ring  source  or  axial  force.  As  a  check  on 
the  program  consider  the  case  of  a  spherical  steel  shell  which 
is  excited  by  a  point  axial  force.  The  following  material  and 
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geometric  constants  in  SI  units  were  used  for  the  calculations: 

Young's  modulus  19.5x10^® 

Poisson's  ratio  0.29 
Density  7700 
Thickness  0.01 
Radius  1.0 

Hysteretic  loss-factor  0.01 
Sound  velocity  in  fluid  1500 
Density  of  fluid  1000 
Distance  force  from  shell  0.5 

Numerical  results  of  far  field  pressure,  dB  ref.  1  micropascal  at 
lm,  are  shown  in  Figure  2.  The  shell  was  idealized  by  41  conical 
shell  elements.  At  low  frequencies  there  is  excellent  agreement 
with  results  obtained  from  a  closed-form  solution  Z21:  as  the 
frequency  increases  the  agreement  worsens,  which  is  typical  of 
results  obtained  by  the  finite  element  method.  Agreement  at  high 
frequencies  could,  no  doubt,  be  improved  by  doubling  the  number  of 
elements  used  in  the  idealization,  but  this  would  be  a  formidable 
computational  task.  Compare  the  numerical  results  with  those  of 
the  point  force  alone;  the  presence  of  the  shell  has  strongly 
enhanced  the  radiation  in  the  mid  frequency  regime.  The  only 
differences  worthy  of  note  between  results  at  0=0  (A  plot)  and 
0=180  (B  plot)  are  the  latter  showing  more  pronounced  radiation 
at  the  340Hz  resonance,  and  reduced  radiation  above  400Hz  because 
the  observation  point  lies  in  the  shadow  zone. 

As  the  source  approaches  the  shell  numerical  results 
(not  shown)  become  more  and  more  inaccurate.  This  is  because  a 
linear  variation  in  pressure  over  an  element,  assumed  in  the 
Helmholtz  integral  equation  discretization  C13,  is  not  suitable 
for  representing  a  surface  pressure  whose  spatial  variation 
becomes  increasingly  like  a  delta  function.  Numerical  checks, 
such  as  doubling  the  number  of  elements,  are  thus  essential. 

As  an  extra  demonstration  of  the  capabilities  of  the 
Fortran  program.  Figure  3  shows  far  field  sound  radiation  when 
the  elastic  scatterer  is  a  prolate  spheroidal  shell  of  length  6m 
and  width  2m,  the  ring  of  radius  lm  being  located  0.5m  aft.  Other 
constants  are  the  same  as  for  the  spherical  shell.  The  point  of 
observation  is  forward  of  the  shell,  viz.  0=180.  In  Figure  3A 
the  excitation  is  a  ring  force;  resonant  enhancement  is  absent, 
but  there  are  frequency  regions  where  the  combined  radiation  of 
the  ring  and  shell  cancel.  In  Figure  3B  the  excitation  is  a  ring 
source;  again  sharp  resonant  enhancement  is  absent,  but  there  is 
a  region  centred  on  250Hz  where  broad-band  enhancement  of  up  to 
5dB  occurs. 

The  sharper  resonant  enhancement  of  the  spherical  shell 
may  be  due  in  part  to  its  modes  having  a  smaller  radiation 
efficiency  than  the  raodeB  of  the  prolate  spheroidal  shell.  Thus 
the  spherical  shell  response  is  not  limited  so  much  by  radiation 
damping.  According  to  Junger  &  Feit  C3D  inefficiently  radiating 
resonances  may  in  certain  circumstances  radiate  more  acoustic 
power  than  their  efficiently  radiating  counterparts. 


4.  CONCLUSIONS 


The  mathematics  of  the  free-space  acoustic  pressure  of  a 
ring  source  and  a  ring  force  have  been  presented  in  a  form  which 
facilitates  computations  of  far  field  sound  radiation  in  the 
presence  of  an  axisymmetric  shell  structure.  Numerical  results 
for  specific  spherical  and  prolate  spheroidal  shells  show  that 
the  former  amplifies  the  radiation,  at  resonance,  much  more  than 
does  the  latter. 

The  programs  can  be  used  to  predict  far  field  sound 
radiation  from  branched  shells  with/without  rib  stiffening,  and 
this  may  be  the  subject  of  a  future  investigation. 


J.H.  James  (PSO) 
November  1988 
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FIG.  2  Far  field  radiation.  Spherical  shell  excited 
by  an  exterior  unit  point  force  (a)  8=0.0  (b)  0=380 
x  x  finite  element  and  Flelrcholtz  integral  equation; 
-  — —  ’exact’; - no  shell. 
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FIG.  3  Far  field  radiation  at  9=180.  Prolate  3  :  1 
spheroidal  shell  excited  by  (a)  exterior  ring  force 
(b)  exterior  ring  source,  x  x  finite  element  and 
Helmholtz  integral  equation;  - no  shell. 
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